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Theoretical  and  Quantum  Chemistry  of  Boron  Compounds 
and  Related  Molecules 


The  theoretical  boron  research  project  under  Contract  Nonr  3471(00) 
for  the  Chemistry  Branch  of  the  Office  of  Naval  Research  has  performed 
research  in  the  development  of  methods  to  describe  the  valence  structure, 
properties  and  behavior  of  boron  compounds  and  related  molecules. 

This  project  has  carried  out  investigations  in: 

I.  Theoretical  Study  of  Charge  Transfer  Behavior  of  Boron- 
Containing  Lewis  Acid-Lewis  Base  Addition  Compounds, 
especially  those  of  the  type  ZB^H^Z. 

II.  Rigorous  Molecular  Orbital  Calculations  of  Boron  Compounds. 

a.  Pentaborane 

b.  Calculation  of  Basic  Parameters  Necessary  for  Quantum 
Chemical  Calculations  on  Boron-Containing  Molecules: 
Modified  (pp  I  pp)  electron  repulsion  integrals  and 
corresponding  boron  screening  constants. 

c.  LCAO-MO-SCF  (Self-Consistent  Field)  Calculations  for 
Non-Saturated  Boron-Containing  Molecules. 

III.  Ionization  Potentials  and  Mass  Spectral  Fragmentation  Patterns 
of  Boron  Compounds;  also  $  v  values  (the  effects  of  sub- 
stitution  on  ionization  potentials)  for  cyclic  radicals  and  for 


O  and  S  compounds 


IV.  Theoretical  Study  of  Aluminum  Coordination  Numbers  in 
Various  Compounds. 

V.  Theoretical  Aspects  of  Bonding  in  Xenon- Fluorine  Compounds. 

Research  on  these  and  other  topics  is  continuing  under  support  from 
the  Air  Force  Office  of  Scientific  Research,  Directorate  of  Engineering 
Sciences,  Propulsion  Division  under  Contract  AF49(638)- 1220. 
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I.  Theoretical  Study  of  Charge  Transfer  Behavior  of  Boron-Containing 
Lewis  Acid-Lewis  Base  Addition  Compounds,  especially  those  of  the 
type  ZB1qH12Z. 

References:  (A-TR-1)  (S-2)  (S-5) 

A  new  class  of  Lewis  acids  and  Lewis  based  -  delocalized  electron- 
deficient  molecules,  such  as  the  higher  borcn  hydrides  -  with  properties 
formally  analogous  to  those  of  delocalized  77-electron  systems  was  de¬ 
fined.  On  this  basis,  the  behavior  of  ZB^H^Z  and  ZB^H^  complexes 
(where  Z  is  an  electron  pair  donor)  could  be  described  in  terms  of  charge 
transfer  complex  theory. 

Treatment  of  the  behavior  of  ZB  -H  Z  compounds  as  a  purely 

XU  i  b 

theoretical  problem  indicated  the  mechanism  by  which  these  complexes 
must  be  formed.  Even  more  important,  derivation  of  the  correct  quantum 
mechanical  description  of  the  valence  structure  and  spectra  of  these 
ZBjqH^Z  compounds  and  subsequent  MO-LCAO  calculations  led  unequivo¬ 
cally  to  the  specific  properties  a  Z  molecule  should  have  in  order  to  be 
most  effective  in  forming  ZB^H^Z  compounds. 

In  PyB^H^Py,  pyridines  act  as  n-donors  to  and  then 

one  gets  reverse  electron  donation  from  B  H  to  the  pyridine  7Tskele- 

XU  lb 

tons.  The  long  wavelength  transitions  in  PyB^H^Py  complexes  were 
shown  to  be  dependent  on  the  energies  of  the  lowest  unfilled  molecular 
orbitals  (LUMO)  of  the  substituted  pyridines  as  calculated  by  MO-LCAO 
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methods  (hence  the  transitions  are  confirmed  as  ■ — *Py).  From 

the  differences  in  the  wavelengths  of  analogous  bands,  the  highest  occupied 
molecular  orbital  (HOMO)  of  B^qH^  would  be  estimated  to  be  0.83  ev 
above  that  of  BQH  .  For  effective  ZB..H  Z  formation  it  is  confirmed 
that  Z  must  not  only  be  a  good  n-donor  but  should  also  be  a  good  7T-( or  d-) 
acceptor. 

Thus  a  purely  theoretical  study  proved  able  to  indicate  to  experi¬ 
mentalists  the  most  judicious  choice  of  which  Z  molecules  should  be  used 
to  form  ZB^H^Z  compounds  and  why  these  particular  Z  molecules 
should  be  the  most  effective.  An  understanding  of  the  theory  of  the  bonding 
in  these  compounds  also  makes  it  possible  to  alter  the  Z  molecule 
correctly  in  order  to  produce  desired  changes  in  behavior 

A  preliminary  report  of  the  application  of  charge  transfer  complex 

(S—  2) 

theory  to  delocalized  electron-deficient  boron  hydrides  was  also  made  '  '. 

II.  Rigorous  Molecular  Orbital  Calculations  of  Boron  Compounds. 

a.  Pentaborane. 

Pentaborane,  B^H^,  was  chosen  as  the  prototype  boron  hy¬ 
dride  molecule  on  which  to  try  several  different  rigorous  and  semi-rigorous 
molecular  orbital  calculational  procedures.  One  of  the  main  purposes  of 
such  a  study  is  to  determine  what  simplifying  approximations  can  be  made 
which  will  most  closely  match  the  results  of  a  rigorous  calculation  and  so 
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will  make  calculations  on  the  higher  boron  hydrides  more  feasible. 

Symmetry  orbitals  for  pentaborane,  B^H^,  were  formulated  for 
several  different  sets  of  basis  orbitals,  both  hybridized  and  unhybridized. 
The  elements  of  the  29x29  Hamiltonian  matrices  have  been  expanded  in 
terms  of  individual  atomic  orbitals.  Integrals  are  being  evaluated  two 
ways:  First,  have  been  estimated  semi-empirically;  Second,  are  being 
calculated  rigorously  by  separating  each  matrix  element  into  kinetic  and 
potential  energy  and  electron  repulsion  terms.  This  has  necessitated  our 
having  to  initiate  setting  up  computer  programs  for  the  integrals  involved 
(specifically  the  many-center  integrals).  A  rigorous  core  approximation 
for  integrals  developed  previously  by  Professor  Daudel  is  being  in¬ 
corporated  into  the  pentaborane  calculations. 

b.  Calculation  of  Basic  Parameters  Necessary  for  Quantum 
Chemical  Calculations  on  Bor  on- Containing  Molecules: 

Modified  (pp  |  pp)  electron  repulsion  integrals  and  corres¬ 
ponding  boron  screening  constants. 

References:  (A-TR-7)  (A-TR-8) 

One  major  improvement  in  MO  calculations  of  boron  compounds 
would  be  the  choice  of  a  more  correct  boron  screening  constant,  for  use 
with  Slater  orbitals,  which  would  embody  electron  correlation  effects. 

Modified  one -center  electron  repulsion  integrals  (ppf  pp) 
for  B  and  B"  which  take  into  account  electron  correlation  effects  have 
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been  evaluated  by  expansion  of  the  energies  of  the  boron  atom  and  its 
positive  and  negative  ions  as  a  function  of  spectroscopic  term  values. 

The  values  for  (pp|  pp)  calculated  from  two  different  pairs  of 
valence  states  7.094  ev  and  6.976  ev  are  very  close.  It  is  demon¬ 
strated  that  for  boron  in  certain  valence  states  the  Pariser-Parr  approxi¬ 
mation  (pp|  pp)  =  I-A  and  the  more  rigorous  approximation  of  Julg 
(ppf  pp)  =  I-A-£(I+A)  do  not  hold,  because  when  the  electron  is  put  into 
a  formerly  empty  orbital,  an  exchange  term  arises  in  addition  to  the 
Coulomb  term. 

For  the  calculation  of  (pp\pp)  -  it  was  necessary  to  evaluate 
the  electron  affinities  of  B  to  form  B~  in  various  valence  states.  This 
was  done  by  adapting  an  existing  rigorous  extrapolation  procedure  to 
valence  state  calculations. 

(pp|pp)g-  was  calculated  by  expansion  of  the  energies  of 
B,  B  and  B  as  a  function  of  spectroscopic  term  values.  The  value 
for  (ppjpp)^-  calculated  in  this  manner,  4.737  ev,  and  the  value  cal¬ 
culated  using  the  Pariser-Parr  approximation,  (ppjpp)  =  I-A  *  4.015  ev, 
are  quite  close.  This  is  especially  gratifying  in  view  of  the  fact  that  the 
two  calculations  are  independent  and  that  each  calculation  depends  on 
extrapolated  values  for  the  electron  affinities  of  B"  in  completely  different 


valence  states 
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In  addition  a  new  general  method  for  calculating  (pp  pp) 
^including  correlation  effects]  for  valence  states  of  any  arbitrary 
atom  was  derived  and  justified  by  a  density  matrix  formalism  for  one 
and  two  electrons  outside  a  "core".  This  general  method  of  evaluating 
(pp  pp)  is  expected  to  be  of  considerable  use  in  all  types  of  quantum 
chemistry  calculations. 

c.  LCAO-MO-SCF  (self-consistent  field)  Calculations  for 
Non-saturated  Boron- Containing  Molecules. 

References:  (A-TR-5)  (Db-3)  (S-7) 


An  LCAO-MO-SCF  (self-consistent  field)  computer  program 

for  non-saturated  systems  has  been  written  and  checked  out  for  7094  IBM 

computer.  The  problem  of  correct  choice  of  input  parameters  for  boron 

(A-TR- 5) 

compounds  was  investigated  '  ,  Calculations  are  now  being  per¬ 

formed  for  boron-containing  molecules  to  obtain  resonance  energies, 
amount  of  partial  double-bond  character,  energy  levels,  charge  dis¬ 
tributions  and  hence  reactivities. 

A  comparative  evaluation  was  made  of  calculations  employ¬ 
ing  only  parameters^  and  [corresponding  to  Coulomb  integrals  (val¬ 
ence  state  ionization  potentials)  and  resonance  integrals]  ,  with  calculations 
including  electron  repulsion  integrals  -  both  with  and  without  iteration  to 
self-consistency.  Boron  ring  compounds  were  used  as  a  test  case  and  ’ 
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several  quite  different  conceptual  models  for  compounds  such  as  borazoles 
were  tried.  These  calculations  were  performed  to  ascertain  how  the 
various  theoretically  calculated  properties  such  as  energy  levels,  charge 
distributions,  localization  energies,  spectra,  etc.  changed  as  a  function 
of  the  calculational  method  employed.  This  problem  was  a  continuation  of 
the  overall  general  problem  which  Dr.  Kaufman  carried  out  at  the  Centre 

i  t 

de  Mecanique  Ondulatoire  Appliquee,  Paris,  France^on  comparison  of 
calculations  by  the  three  techniques  for  all  types  of  molecules.  In  order 
to  be  able  to  evaluate  different  calculational  techniques  for  boron  molecules, 
it  was  felt  necessary  to  investigate  this  problem  for  molecules  for  which 
one  had  much  more  calculational  experience.  Consequently,  during 
May  -  August  1962  these  comparative  calculations  were  carried  out  first 
for  alternant  hydrocarbons,  second  for  non-alternant  hydrocarbons,  then 
for  molecules  containing  heteroatoms.  One  wished  to  see  if  the  geometry 
of  the  molecule  or  the  fact  it  contained  heteroatoms  was  the  dominant  factor 
governing  its  behavior.  A  most  important  finding  was  that  for  sites  of 
chemical  reactivities  there  is  no  consistent  relation  between  charge  dis¬ 
tribution  of  the  parent  molecule  and  localization  energies  of  the  inter¬ 
mediate  complex  (using  several  different  plausible  models  for  the  inter¬ 
mediate  complexes).  This  is  a  point  to  be  carefully  considered  for  re¬ 
activities  of  boron  hydrides,  carboranes,  etc.  where  charge  distributions 

/e  7\ 

are  often  used  as  a  criterion  of  reactivity  . 
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III.  Ionization  Potentials  and  Mass  Spectral  Fragmentation  Patterns 
of  Boron  Compounds;  Also^^  Values  (the  effects  of  substitution 
on  ionization  potentials)  for  cyclic  radicals  and  for  O  and  S  compounds, 

a.  Calculations  of  Ionization  Potentials  and  Mass  Spectral 
Fragmentation  Patterns  of  Boron  Compounds. 

References:  (A-TR-2)  (A-TR-4)  (B-l)  (B-5)  (S-4) 

1)  Mass  spectral  fragmentation  patterns  of  boron  compounds. 

The  monoisotopic  fragmentation  pattern  of  a  molecule  is  very 
important  because  it  provides  great  insight  into  the  mechanism  of  de¬ 
composition  of  the  molecule  under  external  influences  and  is  directly 
related  to  its  thermal  and  radiation  stability.  The  problem  of  calculating 

monoisotopic  fragmentation  patterns  for  boron  compounds  containing  the 
10  11 

natural  B  -B  abundances  was  very  tedious  due  to  the  thousands  of 

separate  terms  which  had  to  be  calculated. 

A  program  for  the  calculation  of  monoisotopic  fragmentation 

patterns  for  boron-containing  molecules  from  mass  spectral  raw  data 

was  written  and  checked  out  with  the  generous  cooperation  of  the  Martin 

Company  Computing  Center.  With  this  program,  it  is  possible  to  calculate 

the  monoisotopic  fragmentation  patterns  of  boron  compounds  labeled  with 

both  B^  and  B**  as  well  as  with  two  isotopes  of  another  element  such  as 

H  or  D.  The  computer  program  is  completely  general  and  is  applicable  to 

any  X  Y  molecule  with  two  different  sets  of  isotopic  labels.  The 
m  n 

availability  of  a  computer  program  to  calculate  monoisotopic  fragmentation 
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patterns  should  make  possible  a  much  wider  use  of  mass  spectrometry 
as  a  means  of  following  the  reactions  and  decomposition  of  boron  compounds 
than  has  been  previously  possible. 

Final  calculations  of  fragmentation  patterns  of  various  isotopically 
labeled  decaboranes  prove  that  the  commonly  used  method  of  estimating 
%  in  a  "normal"  sample  from  the  ratio  of  m/e  10  to  m/e  11  is  NOT 
valid  and  an  alternate  method  must  be  chosen.  This  problem  has  far-reach¬ 
ing  consequences  in  the  boron  field,  especially  when  the  content  of  a 

"normal"  sample  is  counted  on  as  a  neutron  shielding  or  neutron  absorbing 
agent.  The  error  in  estimating  the  %  B10  from  the  ratio  m/e  10  to  m/e  11 
can  be  as  much  as  40%  of  the  total  B1^  content. 

Final  calculations  of  the  monoisotopic  fragmentation  patterns  of 
BgHg  and  B^D^  provide  an  insight  into  the  method  of  decomposition  of 
pentaboranes  which  should  bear  a  relationship  to  their  behavior  in  electric 
discharge  experiments  such  as  are  used  to  prepare  and  car- 

borane  s , 


2)  Calculation  of  ionization  potentials  of  boron  compounds 
from  mass  spectral  appearance  potential  data. 

Calculations  of  the  ionization  potentials  of  B_HQ,  B  D 
and  the  fragments  formed  from  these  molecules  were  completed  during 
this  past  year.  Comparison  of  the  results  indicated  a  very  unusual  re¬ 
verse  H-D  isotopic  effect  on  the  ionization  potentials  which  in  turn  seems 
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to  have  pertinent  implications  with  respect  to  the  valence  structure  of 
pentaborane  and  its  fragments. 

Data  from  an  earlier  appearance  potential  study  of  decaborane 
was  re-examined  to  gain  more  insight  into  the  behavior  of  the  Bjo**12 
fragment.  This  fragment  seems  to  be  quite  closely  related  to 

the  Port*on  ZB^H^Z  Edition  compounds. 

b.  The  Effect  of  Substitution  on  Ionization  Potentials;^  Values. 

K 

References:  (A-TR-3)  (A-TR-9)  (B-2)  (B-3)  (Da-2)  (S-l) 

It  had  been  previously  shown  by  Dr.  Kaufman  that  the  effect  of 

substituent  groups  on  ionization  potentials  of  alkylamines  (Y^N),  alkyl  free 

radicals  (Y^C)  and  boryl  radicals  (Y^B)  were  almost  identical.  From 

measured  photoionization  potentials  of  amines  a  new  set  of  constants.  S 

values,  which  quantitatively  reflect  changes  in  ionization  potentials  were 

derived.  That  thiso^  method  and  the  concepts  on  which  it  is  based  have 

a  general  validity  and  scope  far  beyond  their  use  in  predicting  ionization 

potentials  has  recently  been  confirmed  in  an  article  by  W.  F.  Libby,  J.  Chem. 

Phys.,  35,  1714  (Nov,  1961)  on  radiation  induced  polymerization,  in  which 

our  work  on  &  values  is  cited.  S  v  values  have  also  assumed  importance 

recently  because  of  the  use  of  eT  values  as  a  measure  of  the  magnitude 

of  penetration  integrals  of  neutral  substituent  groups  in  quantum  chemical 
(Da-2)  (S-l) 


calculations. 
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1)  Estimation  of  the  ionization  potentials  of  cycloalkyl 
radicals  and  cyclic  amines  by 6 ^  values; 

( f  „  values  for  alcohols,  ethers  and  sulfides. 

Research  performed  on  this  contract  has  shown  that  cT  values 
obtained  from  linear  and  branched  alkyl  substituents  proved  to  be  extendable 
to  cyclic  substituents  in  which  the  carbon  from  whence  the  unpaired  electron 
was  being  withdrawn  was  itself  part  of  the  ring.  This  fact  would  not  have 
been  apparent  as  simply  a  logical  extension  of  our  earlier  work  on.^> 

K 

value  s . 

Effects  of  substitution  on  ionization  potentials  were  proved  not  to 
be  the  same  for  the  series  of  alcohols  and  ethers  as  for  the  amines,  a 
specific  set  of  values  for  use  with  oxygen  compounds  was  derived  from 

photoionization  measurements.  The  purpose  of  these  values  is  twofold: 
first,  to  enable  one  to  estimate  ionization  potentials  of  other  oxygen- 

containing  molecules  and  second,  to  permit  one  to  differentiate  effects 

/  /  • 

of  penetration  integrals  on  different  cores,  -C  and  -N  :  on  one  hand, 
and  -O-  on  the  other  hand.  An  analogous  set  of  values  for  S 
compounds  proved  to  be  different  from  that  for  O  compounds  and  verifies 
the  fact  that  the  effect  of  substituent  groups  is  different  on  first  and  second 
row  cores.  It  is  concluded  that  one  should  use  different  magnitudes  for 
penetration  integrals  for  substituents  on  first  and  second  row  elements. 

The  observed  differences  in  substituent  group  effects  depending  on  whether 
an  O  (or  S)  radical  or  molecule  is  being  ionized  imply  that  one  should 


perhaps  use  different  evaluations  of  penetration  integrals  depending  upon 
whether  one  or  two  electrons  from  O  (or  S)  are  contributed  to  a  delocali- 
zed  ,y..em.<A-TR-3)  <B-2>  ‘B-3> 


2)  Extension  of  the  density  matrix  treatment  to  a  two- 
body  model  for  electrons  outside  a  "core". 

Extension  has  been  made  to  a  two-body  model  for  a  pair  of 

electrons  outside  a  "core"  of  the  treatment  by  density  matrix  formalism 

of  a  one-body  model  for  an  electron  outside  a  "core".  S  values 

K 

[differences  in  ionization  potentials  between  the  parent  compounds] 

for  the  7Telectron  of  alkyl  radicals  can  be  developed  exactly  as  a  one -body 

model  problem  in  terms  of  density  matrices.  The  derivation  of  the  two-body 

model  permits  one  quantum  mechanically  to  describe  the  same  &  effects 

K 

on  the  ionization  potentials  of  an  electron  from  the  lone  pair  of  substituted 
amines  by  an  analogous  treatment  using  the  same  perturbation  Hamiltonian 
as  for  substituted  alkyl  radicals  with  two  electron  functions  for  the  lone 


pair  of  the  nitrogen. 


(A-TR-9) 


IV.  Theoretical  Study  of  Bonding  in  Aluminum  Compounds. 


A  logical  extension  of  our  interest  in  boron  compounds  led  us  to 
examine  analogous  properties  of  aluminum  compounds.  The  theory  of 
charge  transfer  complexes  was  applied  to  the  aluminum  compounds  and  on 
this  basis  a  great  deal  of  regularity  becomes  apparent  in  the  properties  of 
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alumiaum  compounds  with  similar  coordination  numbers:  tetra-,  penta- 
or  hexacoordinate.  Infrared  and  NMR  data  on  these  compounds  were 
examined  and  criteria  for  the  determination  of  the  aluminum  coordination 
numbers  in  various  compounds  were  formulated. 

Definition  used  for  electron  donors  in  Mulliken  charge -transfer 
theory  are  invoked  in  describing  coordination  to  aluminum.  Z,  a 
neutral  electron  pair  donor,  and  Y  ,  an  ion  electron  pair  donor,  are 
considered  as  equivalent. 

Aluminum  NMR  absorptions  should  exhibit  the  same  shift  to 

higher  field  with  increased  A1  cotirdinacy  as  do  boron  compounds.  In 

addition,  the  magnitude  of  the  coupling  constant  between  hydrogen  and 

aluminum  should  depend  significantly  on  the  percentage  of  s -character  of 

the  bonding  hybrid  of  the  directly  bonded  nucleus.  NMR  line  widths  of  A1 

compounds  have  also  been  considered  indicative  of  the  cubic  (octahedral 

27 

or  tetrahedral)  or  noncubic  symmetry  around  the  A1  nucleus. 

An  Al-H  stretch  observed  in  the  infrared  has  been  reported  to 
be  indicative  of  whether  an  A1  atom  is  tetra  -  (5.  6k)  or  pentacodrdinate  - 
(5.8  -  5.9K).  This  shift  is  easily  understood,  since  the  force  constant  of 
a  bond  to  a  tetracodrdinate  A1  must’ be  stronger  than  that  of  a  bond  to  a 
pentacodrdinate  Al,  and'K^ykT 

All  of  these  criteria  should  provide  an  indication  of  the  coordination 
number  of  Al  in  a  specific  compound,  and  several  examples  are  discussed. 
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V.  Theoretical  Aspects  of  Bonding  in  Xenon-Fluorine  Compounds 
(A-TR-10). 


Theoretical  concepts  of  bonding  in  xenon-fluorine  compounds  were 
developed  to  permit  one  to  generalize  and  to  predict  the  existence  or 
non-existence  of  other  closely  related  rare-gas  compounds. 

The  combination  of  rehybridization  of  the  xenon  atom  to  utilize 
its  d  orbitals  in  valence -shell  bonding  and  back-donation  of  TFelectrons 
from  fluorine  is  the  factor  which  contributes  to  the  formation  and  stability 
of  these  Xe-F  compounds. 

The  correct  geometry  for  XeF^  was  predicted  from 
these  postulates,  as  well  as  the  identification  of  the  lower  xenon- 
fluoride  as  XeF^  and  its  correct  geometry,  and  the  existence  of  both 


Kr-F  and  Rn-F  compounds 
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